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This paper describes the numerical study on film cooling effectiveness and aerodynamic loss due to coolant and 


main stream mixing for a turbine guide vane. The effects of blowing ratio, mainstream Mach number, surface 


curvature on the cooling effectiveness and mixing loss were studied and discussed. The numerical results show 


that the distributions of film cooling effectiveness on the suction surface and pressure surface at the same blowing 


ratio (BR) are different due to local surface curvature and pressure gradient. The aerodynamic loss features for 


film holes on the pressure surface are also different from film holes on the suction surface. 
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Introduction 


Film cooling is one of the most effective technologies 
to protect turbine vanes and blades from the hot gas. Pre- 
vious film cooling studies focused on the influences of 
geometric and flow parameters on the adiabatic cooling 
effectiveness [1,2]. The effects of blowing ratio, density 
ratio, turbulence intensity, hole incline/compound angle, 
length-diameter ratio and hole shapes on the film cooling 
effectiveness were well documented [3,4]. Film cooling 
reduces the turbine metal temperature. Meanwhile, the 


mixing of coolant and hot gas leads to kinetic energy loss. 


The additional aerodynamic loss (or gain) due to the 
coolant and main stream mixing has a strong impact on 
the turbine thermal efficiency, especially for modern gas 
turbines in which full coverage film cooling is commonly 
applied. In the study of film cooling characteristics, the 
effects of film coolant ejection on turbine aerodynamic 
loss should also be taken into account. 

Haller et al. [5] studied the effects of coolant ejection 
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on the loss of vane cascade. The experiments showed that 
the aerodynamic loss due to coolant and main stream 
mixing is related to the locations of film holes, blowing 
ratios and exit Mach number, and the mixing loss is in- 
creased with the increase of blowing ratio. Day et al. [6] 
performed experiments to study the effects of coolant 
ejection on the aerodynamic loss of a turbine guide vane. 
The experimental results showed that the cascade aero- 
dynamic loss depends on the coolant ejection location, 
number of film hole rows and film hole shape. The mix- 
ing loss due to fan-shaped holes is generally higher than 
that of round holes. The mixing loss due to the coolant 
ejection from the film hole on the suction surface is 
higher than that on the pressure surface. Gomes et al. [7] 
studied the effects of coolant ejection on flow separation 
on suction surface of a high loading turbine vane. It is 
found that the flow separation appearing at low Reynolds 
number can be restrained by the coolant ejection and the 
total aerodynamic loss can be reduced. Friedrichs and 
Hodson [8] experimentally studied the effects of endwall 
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Nomenclature 


BR Blowing ratio, - 


D Diameter of cooling hole, mm 

H Pitch of cooling holes, mm 

I Momentum ratio, - 

K Kinetic energy ratio, - 

Laxial Axial chord length of the vane, mm 
Le Chord length of the vane, mm 

m Mass flow rate, kg/s 

Mex Mach number at cascade exit, - 

P Static pressure, Pa 

P; Total pressure, Pa 

S Curve length to the cooling hole center, mm 
T Total temperature, K 

4 Velocity, m/s 


film cooling on the main flow. They found that coolant 
ejection can change the secondary flow structure in the 
main flow passage and even reduce the loss of secondary 
flow. The aerodynamic loss is enlarged with the increase 
of blowing ratio for the endwall film cooling. Colban and 
Thole [9] studied the effects of endwall film hole shape 
on the aerodynamic loss. They found that film cooling 
with fan-shaped hole results in smaller loss than with 
round film hole. Rehder [10] studied the effects of coo- 
lant trailing-edge ejection on aerodynamic loss. It was 
found that the aerodynamic loss is small while the coo- 
lant is ejected from the slot at trailing edge center, and it 
is large when the coolant is ejected from the slot at pres- 
sure surface. Schobeiri and Pappu [11] investigated the 
correlations between the type of trailing edge slot and 
coolant to main stream momentum ratio, velocity ratio 
and trailing edge thickness. There are optimal momentum 
ratio, velocity ratio and trailing edge that lead to mini- 
mum aerodynamic loss. 

Most of previous studies focused on either film cool- 
ing effectiveness or aerodynamic loss. In this paper, nu- 
merical simulations were carried out to study both the 
film cooling effectiveness and the aerodynamic loss due 
to coolant and main stream mixing for a turbine guide 
vane. The effects of blowing ratio, main stream Mach 
number, surface curvature on the cooling effectiveness 
and mixing loss are presented and discussed. After this 
introduction, numerical method used for the simulation, 
validation of the method and main results are presented. 


Numerical Method and Validation 


Ansys CFX was employed to simulate the main-flow 
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a Injection angle, ° 

Naw Adiabatic cooling effectiveness, - 
P Density, kg/m? 

Čp Total pressure loss coefficient, - 


Total pressure loss coefficient generated inside 


CPhole the cooling hole, - 

) Total pressure loss coefficient generated by 
CPmix mixing and inside boundary layer, - 

l Total pressure loss coefficient for film cool- 
ČPiotal in g, - 

Subscripts 

c Coolant 

m Mainstream 

cin Cooling hole inlet 

cout Cooling hole outlet 

w Wall 

f Computational domain outlet 


and coolant flow in a turbine guide vane with film cool- 
ing. The solver is based on finite volume method, and is 
second order accurate. Fully implicit method is used for 
the discretization of Reynolds averaged Navier-Stokes 
equations. 

SST turbulence model with scaled wall function was 
employed for the present simulations, which can 
smoothly change from the solution of RANS equation to 
wall function near the wall. SST model provided highly 
accurate prediction of flow separation under adverse 
pressure gradients, which are common in the boundary 
layer on the suction surface near the trailing edge. As 
film cooling was introduced, several studies [12-14] 
showed that the SST turbulence model is suitable for the 
film cooling calculations. 

Convergence to a steady state is reached when the 
overall residuals of the primary variables are less than 1x 
10% and the variables at monitored points are kept un- 
change for enough iteration steps. 

A validation simulation was performed for the case 
reported in [15], and then compared with the experiment 
result of Ito et al. [16], as shown in Figure 1. For the 
pressure surface film cooling, calculated loss profiles 
agree well with the experiment results. For the suction 
surface film cooling, calculated loss profiles at BR=1.0 
agree well with the experiments at BR=0.92, except 
small discrepancies at y/p=0 and y/p=0.2. While at BR= 
2.0, the calculated loss profiles were slightly higher than 
the experiments in the middle of the trailing edge wake. 
Both calculated loss profiles show the same variation 
trends as the experimental ones. The comparisons of nu- 
merical results and experiment results show that the nu- 
merical method applied is feasible and reliable. 
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Calculation Setup and Boundary Conditions 
Geometry and grid 


The turbine guide vane with a row of film cooling 
holes on both suction and pressure surface is shown in 
Figure 2. The overall chord (L.) of the vane is 180 mm 
and the axial chord (Loyiai) is 94.1 mm. The cylindrical 
cooling hole is located at 56 percent (suction surface) and 
22 percent (pressure surface) of the axial chord down- 
stream of the airfoil leading edge, respectively. Diameter 
of the cooling holes is D=2 mm, with an injection angle 
of a=23°. The pitch of the cooling holes (H) is 5D and 
the vane-to-vane pitch of the cascade is 150 mm. 

The computational domain includes the cascade flow 
passage, the cooling hole and the coolant plenum. The 
height of one hole pitch (H=5D) at middle span of the 
cascade was chosen to decrease the grid scale and avoid 
the influence of secondary flow at the endwall. Periodic 
conditions were applied between the hub and shroud sur- 
face to model an infinite long vane. Coolant plenum re- 
gion was extend to 10D in both direction to simulate the 
full size coolant plenum. 
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Fig. 1 Comparisons of total pressure loss coefficient profiles 
between computed and measured results 
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Fig. 2 Turbine cascade with film cooling holes and computa- 
tional domain 
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Fig.3 Details of the grid near film hole: (a) Symmetry plane of 
the cooling hole; (b) Surface mesh of plenum and vane 


Three different structured grids were generated re- 
spectively for vane cascade without film holes (NO 
HOLE), vane cascade with cooling holes on the suction 
surface (SSH) and vane cascade with cooling holes on 
the pressure surface (PSH) by using ICEM CFD. Double 
O-block grid was implemented in cooling holes to im- 
prove the local grid quality, as shown in Figure 3. The 
first layer height near the wall was 0.01 mm and the 
growth ratio was 1.2. Thus, Y° was less than 3 for all 
grids near solid walls. 

As shown in Figure 4, the cascade flow domain was 
divided into 31 blocks so that the local grid quality can 
be controlled separately. Different scale of the grid can be 
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generated by change node numbers in the blocks. The 
grids of leading edge block 114 and the trailing edge 
blocks 74 and 142 were refined to improve the prediction 
accuracy in these areas. The blocks downstream of the 
cooling hole was refined as well, to make sure that the 
mixing between coolant and main flow can be captured 
accurately. 

Grids for both SSH and PSH contain approximately 
5.7 M nodes. A SSH grid of 6.3M nodes was generated 


by refining the special blocks to check grid independence. 


No visible change was obtained for both cooling effec- 
tiveness and downstream profile loss, as shown in Figure 
5. So it can be judged that the present simulation was 
grid independent. 
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Fig. 4 Mesh block details for cascade flow domain 


Boundary Conditions 


As mentioned above, periodic conditions were applied 
to the hub-shroud interfaces to model an infinite linear 
cascade, and the vane-to-vane interface was also set as 
periodic boundary. On the wall surfaces of plenum, cool- 
ing hole and vane, adiabatic no slip wall conditions were 
applied. At the outlet of the computational domain, static 
pressure was specified as P=1.01 atm. Main stream and 
coolant were both assumed as perfect gas. The total tem- 
perature of main stream and coolant were 500K and 
333K, respectively. Various inlet velocity or total pres- 
sure was applied to different cases to obtain different BR 
and/or Mex. The inlet turbulence intensity was 5% for 
both the main stream and the coolant. 


Test Cases 


There were 15 test cases calculated in the present 
work with various inlet conditions and blowing ratios, as 
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shown in Table 1. Case 0 was a vane cascade without 
cooling hole (NO HOLE), and was a reference case. The 
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Fig.5 Grid independence verification 


Table 1 Test Cases 


Hole Position BR Mex Notation 
Case 0 - - 0.17 NO HOLE 
Case | SS 0.5 0.17 SSH 
Case 2 SS 1.0 0.17 SSH 
Case 3 SS 1.5 0.17 SSH 
Case 4 SS 2.0 0.17 SSH 
Case 5 PS 0.5 0.17 PSH 
Case 6 PS 1.0 0.17 PSH 
Case 7 PS 1.5 0.17 PSH 
Case 8 PS 2.0 0.17 PSH 
Case 9 SS 1.0 0.3 SSH 
Case 10 SS 1.0 0.5 SSH 
Case 11 SS 1.0 0.7 SSH 
Case 12 PS 1.0 0.3 PSH 
Case 13 PS 1.0 0.5 PSH 
Case 14 PS 1.0 0.7 PSH 
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inlet velocity of Case 0 was 20m/s, the corresponding 
cascade exit Mach number was 0.17. This inlet velocity 
was also applied to cases 1 to 8. Cooling hole for cases 1 
to 4 was on the suction surface (SSH) and was on the 
pressure surface for cases 5 to 8 (PSH). Cases 9 to 14 
maintained the blowing ratio as an invariant and various 
exit Mach number were obtained from 0.3 to 0.7 by set- 
ting different inlet total pressure at main stream inlet. 


Parameter Definitions 


Define blow ratio, momentum ratio and kinetic energy 
ratio of the coolant to the main stream as follows, 


Jee (1) 


2 
l= PV (2) 


Klee (3) 


Where, p, and V, were the coolant density and velocity, 
Pm and V,, were the main stream density and velocity. 
Define the adiabatic film cooling effectiveness as, 


=m? Ww (4) 


Where T„ and T, were the inlet temperature of main 
stream and coolant. T„ was the adiabatic temperature on 
the vane surface. 

For the vane without film cooling hole, the total pres- 
sure loss coefficient čp yane was defined as, 
= Fam E P, f 


ane = 5 
SP aie Pn -P; ( ) 


Where P,m was the main stream inlet total pressure, 
and P,, and P; were the outlet total pressure and static 
pressure, respectively. 

For the cooled vane, total pressure loss coefficient 
should be transformed as: 


m, Mn 
: 7 P cin to 7 Fam j Pf 
E = M, +My, mM, +My (6) 
Ptotal P P 
tm * f 


Where P,cin was the total pressure at the inlet of the 
cooling hole. m, and m,, were respectively the mass 


flow rates of the coolant and the main stream. The total 
losses include the loss in cooling hole, the loss due to 
mixing process and the loss in cascade boundary layer. If 
replacing P, cin with the total pressure P, cous at the cooling 
hole outlet, we can obtain the total pressure loss coefficient 
ČP mixing that contained the latter two parts of čp rotal 1.€., 

m m 


c P + m P_-P, 
. : t,cout . . t,m t,f 
E = mM, +My M, +My (7) 
P,mix P P 
tm of 
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Obviously, the loss inside the cooling hole can be ob- 
tained as follows, 
m 
pag een ~ P cout) 
ŠP hole T ŠP total 7 SP mix =—* (8) 
Lg Ep 


tm 
Results and Discussions 


Cooling Effectiveness at Different BR 


To reduce the vane adiabatic temperature was the im- 
portant purpose of film cooling. Cooling Effectiveness 
was an indicator of temperature reduction by using film 
cooling. Figure 6 showed the lateral averaged adiabatic 
film cooling effectiveness on the suction surface and the 
pressure surface at different BR (Cases | to 8). S/D is the 
ratio of curve length from the center of the cooling hole 
outlet to the hole diameter. For film cooling on suction 
surface (SSH), cooling effectiveness at BR=0.5 is the 
highest at almost all district downstream of the cooling 
hole. While for film cooling on pressure surface (PSH), 
the cooling effectiveness at BR=0.5 is the highest near 
the cooling hole outlet and then dropped faster than that 
of SSH. The different cooling features between SSH and 
PSH may be caused by the surface curvature, normal 
pressure gradient and also stream-wise pressure gradient. 

For both SSH and PSH, the same dropping trend of 
cooling effectiveness was found when BR was increased. 
This can be explained as that the high normal momentum 
of the coolant raised the coolant stream apart from the 
vane surface further. With the increase of BR, the cooling 
effectiveness for SSH dropped much faster than for PSH. 
For cases 6 to 8, the cooling effectiveness decreased till 
S/D=35, and then rose up downstream, as a result of the 
coolant reattachment due to the concave curvature of the 
pressure surface. 
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Fig. 6 Lateral averaged adiabatic cooling effectiveness for 
SSH and PSH at different BR (M.,=0.17) 
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Cooling Effectiveness at Different Mex 


The lateral averaged adiabatic cooling effectiveness 
profiles for SSH and PSH at BR=1.0 and different Mex 
were shown in Figure 7. For SSH, the cooling effective- 
ness downstream of the cooling hole increased as the exit 
Mach number varied from 0.17 to 0.5. Compared with 
Me = 0.5, the cooling effectiveness at Mex = 0.7 was 
higher at S/D<20, and then almost identical at S/D>20. 
The cooling effectiveness for SSH is increased from the 
cooling hole exit to the trailing edge. The coolant jet at 
BR=1.0 was slightly raised up from vane suction surface, 
but vane surface could still be cooled by the coolant 
stream. The same profile trend was found for SSH at 
BR=1.0 in Figure 6. 

Reattachment of coolant stream to the pressure surface 
was found for PSH at different M.,. The effects of coo- 
lant reattachment on the film cooling effectiveness are 
obviously increased as Mex is increased. At Mex = 0.5 and 
M.x= 0.7, the cooling effectiveness rose sharply to 0.15 
and 0.23 near the trailing edge, because the high speed 
coolant stream struck onto the pressure surface near the 
vane trailing edge and then spread in the span-wise of the 
vane surface. To illuminate the reattachment of coolant, 
Figure 8 showed the total temperature contours plotted 
on five planes perpendicular to PS at different x/Z,,;q; for 
the case of Mex = 0.7. It can be seen that main stream 
pushed the coolant onto the pressure surface. The core of 
coolant was forced to spread in the span-wise direction. 


Aerodynamic Loss at Different BR 


Total pressure loss coefficient for SSH at different BR 
was shown in Figure 9. All Epnote, €pmix and Čprotai Were 
increased while BR was changed from 0.5 to 2.0. The 
increment of Epnoie Was approximately the same as the 
increment of &p,,;.. As a result, the percentage of čprole IN 
Eprotals Was increased. At higher BR, the loss generated in 
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Fig. 7 Lateral averaged adiabatic cooling effectiveness for 
SSH and PSH at different Mex (BR=1.0) 


the cooling hole became more important and was not 
negligible. 

The flow distributions inside the cooling hole at BR= 
0.5 and BR=2.0 were show in Figure 10. Separation re- 
gion in the cooling hole located near the hole inlet. 
Another separation region out of the cooling hole was 
found near the hole outlet. The coolant stream in the 
cooling hole separated much more severely at the high 
BR than at the low BR. 

The boundary layer downstream of the cooling hole 
was thinned in the middle plane between the two holes as 
BR increased, as shown in Figure 11 (a). In the meri- 
dional plane of the computational domain, different trend 
was found as shown in Figure 11 (b). At BR=0.5 and 
BR=1.0, the boundary layer was thickened by the coolant 
jets downstream the cooling hole. At BR=1.5 the coolant 
rose apart from the suction surface and the boundary 
layer was thinned. At BR=2.0, the boundary layer thick- 
ness increased because of the severe separation near the 
cooling hole outlet. 

Figure 12 showed the total pressure loss coefficient 
profiles at the domain exit at different BR. y/p>0 was for 
the suction surface and y/p<0 was for the pressure sur- 
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Fig. 8 Total temperature contours on five different planes 
perpendicular to pressure surface (BR=1.0, Mex=0.7) 
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Fig. 9 Total pressure loss coefficient for SSH at different BR 
(Ma=0.17) 
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(b) BR=2.0 


Fig. 10 Flow distribution inside the cooling hole for SSH 
(M,x=0.17) 


face. It was obvious that loss coefficient profile for SSH 
was the highest at BR=2.0 and was the lowest at BR=0.5, 
which was coincident with Figure 9. As BR was in- 
creased, the baseline of the loss rose up and the differ- 
ence between the peak and the baseline was decreased. 
As a result, the loss profiles became flatter at high BR. 

The locations of loss peaks for the cooled vane def- 
lected to the suction surface compared with the uncooled 
vane, because the SSH coolant stream contained a nor- 
mal velocity component towards the mainstream, which 
gave the wake a slight offset to the normal direction of 
the vane surface. Then the offset expanded to the exit and 
a visual deflection was found on the loss profiles. 

The same analysis was implemented for PSH. Table 2 
showed the total pressure loss coefficient for PSH at dif- 
ferent BR at M,,=0.17. Epnoie Was increased two orders of 
magnitude, and was still negligible compared with Cp mix, 
which was the main portion in €pjoq;. The reason for the 
less difference between PSH loss and NO HOLE loss 
was discovered by further investigation. Coolant stream 
of PSH stuck in the boundary layer of pressure surface 
and less impact was caused on the mainstream. As a re- 
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Fig. 11 Boundary layer flow for SSH (M,,=0.17) 
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Fig. 12 Total pressure loss coefficient profiles for SSH at 
different BR (M.,=0.17) 


Table 2 Total pressure loss coefficient for PSH 


ČPhole CR mix CProtal 
NO HOLE - - 4.57E-02 
BR=0.5 5.00E-06 4.57E-02 4.57E-02 
BR=1.0 3.08E-05 4.57E-02 4.57E-02 
BR=1.5 9.13E-05 4.57E-02 4.58E-02 
BR=2.0 2.00E-04 4.59E-02 4.60E-02 
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sult, Epioai increased slightly with the increase of BR, 
because pm maintained almost unchanged. Separation 
in the cooling hole and thickness increment in the boun- 
dary layer were observed at the same location and the 
intensity of the separation had the same variation among 
different BR. 

Figure 13 showed the relations between the loss coef- 
ficient and the momentum ratio I or the kinetic energy 
ratio K. For Sproles €pmix, ANd Čprotai , the relation was non- 
linear with the momentum ratio I and was linear with the 
kinetic energy ratio K. So the aerodynamic loss caused 
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Fig. 13 Relations between the loss coefficient and momentum 
ratio or kinetic energy ratio (M.,=0.17) 
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by film cooling was in proportion to the third power of 
the velocity ratio of coolant to main stream. The same trend 
was found for PSH, as shown in Figure 13 (b) and (c). 

An equation can be fitted as, 


3 
E,=A-K+B=ALOe +B (9) 

Pml in 
In which A was the curvature factor and B was the in- 
tercept at BR=0. The values of A and B for SSH and PSH 
were listed in Table 3 and Table 4. This equation can be 
used to predict the total pressure loss coefficient when 

the density ratio and velocity ratio were known. 


Aerodynamic Loss at Different Mex 


Total pressure loss coefficients at different Mex (BR= 
1.0) were shown in Figure 14. Čpmix and €pjora for both 
SSH and PSH shared the same trend, i.e., it is slightly 
decreased with the increase of M,, at Mex<0.4, and then 
increased at M,,>0.4. There was a proper M,, to gain the 
minimum aerodynamic loss. Inner hole loss coefficient 
Epnote for SSH dropped as Mex increased. The values of 
ČPhole for PSH and SSH shown in Figure 14 were listed in 
Table 5. 

Total pressure loss coefficient profiles at different Mex 
was shown in Figure 15. Deflections of the loss peaks 
towards y/p<0 were found for both SSH and PSH at high 


Table 3 Values of factor A in equation 9 


CPhole op mix Čr total 
SSH 7.78E-03 8.61E-03 1.64E-02 
PSH 3.44E-05 3.64E-05 6.65E-05 


Table 4 Values of factor B in equation 9 


CPhole CP mix ČPtotal 
SSH 4.01E-04 4.35E-02 4.39E-02 
PSH 5.15E-06 4.56E-02 4.57E-02 
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Fig. 14 Total pressure loss coefficient for SSH and PSH at 
different Mex (BR=1.0) 


58 


Table 5 Inner hole loss coefficient at different Mex 

M.=0.17 M..=0.30 M..=0.50 M.,=0.70 
SSH 6.14E-03 5.13E-03 4.51E-03 3.98E-03 
PSH 3.08E-05 2.44E-05 2.10E-05 1.77E-05 

0.25- — SSH M, =0.17 

| Fon bee = SSHM,,=0.30 

0.204 / A - - - -SSH M,„=0.50 

| YN |+---- SSH M, =0.70 
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- | 
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0.05 4 
0.00 = 
a eee en eel 
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Fig. 15 Total pressure loss coefficient profiles at different Mex 
(BR=1.0) 


Mex. As Mex increased, main stream speed dominated the 
deflection because of the vane curvature. Different from 
the changes caused by BR, the peak value in the profiles 
increased and the baseline value decreased as Mex in- 
creased, which meant that the losses concentrated in the 
region influenced by trailing edge. The increment of 
main stream speed leads to the increase of friction loss 
along the vane surface and mixing loss between coolant 
stream and main stream. 


Conclusions 
A validated numerical method was used to study cy- 


lindrical hole film cooling effectiveness and aerodynamic 
loss features of a turbine guide vane at different BR and 
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different Mex. The trends of cooling effectiveness varia- 
tion were explained and the loss variation mechanisms 
were elucidated. Specific conclusions are drawn as fol- 
lows: 

(1) For both SSH and PSH, the highest cooling effec- 
tiveness is gained when cooling holes are located at suc- 
tion surface at BR=0.5 (M,,=0.17). Mach number of 
main stream has an impact on the film cooling effective- 
ness. In general, averaged cooling effectiveness is in- 
creased with the increase of Mex. Rapid increase of cool- 
ing effectiveness for PSH at high Mex was noticed. 

(2) Total pressure loss coefficients Epnoie, Epmix, and 
Eprom Were in linear relation with the kinetic energy ratio. 
A linear equation was fitted for both SSH and PSH. Total 
pressure loss coefficients are also slightly impacted by 
Mex, as well as BR (I or K). Loss coefficient profile at 
domain outlet became flatter as BR is increased. 
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